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Abstract. During the last few years, the phase diagram of the large N Gross-Neveu 
model in 1 4-1 dimensions at finite temperature and chemical potential has undergone a 
major revision. Here we present a streamlined account of this development, collecting 
the most important results. Quasi-one-dimensional condensed matter systems like 
conducting polymers provide real physical systems which can be approximately 
described by the Gross-Neveu model and have played some role in establishing its phase 
structure. The kink-antikink phase found at low temperatures is closely related to 
inhomogeneous superconductors in the Larkin-Ovchinnikov-Fulde-Ferrell phase. With 
the complete phase diagram at hand, the Gross-Neveu model can now serve as a firm 
testing ground for new algorithms and theoretical ideas. 



Updating the Gross-Neveu phase diagram 



2 



1. Introduction 

In its original form, the Gross-Neveu (GN) model PP is a relativistic, renormalizable 
quantum field theory (QFT) of N species of self-interacting fermions in 1+1 dimensions 
with Lagrangian 



The bare mass term ~ itlq explicitly breaks the discrete chiral symmetry ip 'j^ip of 
the massless model. The interaction term can be generalized to acquire a continuous 
chiral symmetry as in the Nambu-Jona-Lasinio (NJL) model, but we shall not consider 
this option here. As far as the phase diagram is concerned, the 't Hooft limit N ^ oo, 
g^ ~ is most instructive. The reason is the following. In 1+1 dimensions, classic no- 
go theorems forbid spontaneous breaking of a continuous symmetry at zero temperature 
or of a discrete symmetry at finite temperature [21 El- Although mean-field theory 
often predicts these effects, fluctuations are expected to destroy any long range order. 
Provided that one deflnes the model by letting iV — oo before taking the thermodynamic 
limit, these fluctuations are suppressed and a richer phase structure becomes accessible 
We follow this line of reasoning here. 

Like in a previous pedagogical review article [S] which the present paper is meant 
to update, we work canonically, couching the theory in the language of the relativistic 
Hartree-Fock (HP) approximation. It is formally equivalent to the semi-classical 
functional integral approach but conceptually simpler. Physically, it emphasizes the 
fact that the Dirac sea is an interacting many-fermion system [HIIZI, thereby making it 
somewhat easier to switch back and forth between the relativistic QPT and condensed 
matter physics literature. 

Let us now briefly recall the motivation for studying the Lagrangian (Q). QPT 
models for which one can construct the exact renormalized phase diagram are extremely 
scarce and worth investigating on theoretical grounds, even if they are far from being 
realistic. As a matter of fact, in spite of its simple Lagrangian, the GN model 
shares many non-trivial properties with quantum chromodynamics (QCD), notably 
asymptotic freedom, dimensional transmutation, meson and baryon bound states, chiral 
symmetry breaking in the vacuum as well as its restoration at high temperature and 
density. Perhaps even more surprising and less widely appreciated is the fact that 
GN type models have enjoyed considerable success in describing a variety of quasi-one- 
dimensional condensed matter systems such as the Peierls-Prohlich model, conducting 
polymers like polyacetylene, or inhomogeneous superconductors. 

The title of this paper has the following background: Originally, the kink and 
kink-antikink baryons flrst derived in QPT by Dashen, Hasslacher and Neveu |H] have 
been useful for understanding the role of solitons and polarons in electrical conductivity 
properties of doped polymers [5]. In the following years, a lot of progress was made 
in condensed matter theory towards understanding polaron crystal structures. This in 
turn helped us to construct the full phase diagram of the relativistic QPT. However, in 
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order not to confuse the issues, we shall first discuss the GN model without reference 
to condensed matter physics and comment on the relationship only towards the end. 

Our last remark about the motivation concerns the restriction to 1+1 dimensions. 
At first glance, this looks very unrealistic. However, most of the work discussed below 
has to do with finite chemical potential. In the presence of a Fermi surface, a kind of 
dimensional reduction takes place which has been exploited for instance in high density 
effective theory (HDET) to QCD (see [ID] for a review). It also manifests itself through 
the well-known fact that Cooper pairing in superconductivity occurs for arbitrarily weak 
attraction, just like in one space dimension. Hence, from the physics point of view, the 
restriction to 1+1 dimensions is perhaps a better idea in the presence of fermionic matter 
than naively thought. 

The Lagrangian has two bare parameters, g"^ and mo. In the process of 
regularization and renormalization, all observables can be expressed in terms of two 
physical parameters m and 7. The relation to the bare quantities and the ultra-violet 
(UV) cutoff A is given by the vacuum gap equation 

-— = 7 + ln-, 7:=-—^. 2 

A(yf^ m JS g"^ m 

Whereas the physical fermion mass m in the vacuum merely provides the overall mass 
scale and can be set equal to 1, the second parameter 7 (called "confinement parameter" 
in condensed matter physics) parametrizes different physical theories. It measures the 
amount of explicit chiral symmetry breaking and vanishes in the massless (mo = 0) 
limit. Notice that it can also be expressed in terms of the physical fermion masses at 
bare mass mo and in the chiral limit. 

The main subject of this paper is the full phase diagram of the massive GN model 
as a function of temperature T, chemical potential fi and confinement parameter 7. The 
original works date from 1985 for the massless ^JJ and 1995 for the massive case J2] 
(see also ^2] for earlier, partial results). As has become clear by now, the assumption 
that the scalar condensate (ipip) is spatially homogeneous made in these works is too 
restrictive. It misses important physics closely related to the existence of kink-antikink 
baryons in the GN model. In the present paper we describe how this problem has 
been cured. This is not a small correction, but requires a substantial new effort. A 
homogeneous condensate acts like a mass and reduces the thermodynamics essentially to 
that of a massive, free relativistic Fermi gas, the mass being determined self-consistently. 
By contrast, solving the Dirac equation with a periodic potential and establishing self- 
consistency is a highly non-trivial task, even in 1+1 dimensions. Due to a number of 
lucky circumstances, this is nevertheless possible in the case at hand, to a large extent 
even analytically. 

In our recent papers we had to proceed from special cases to more general ones 
(treating mo = before mo 7^ 0, T = before T 7^ 0) in order to reduce the problem 
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to manageable size. Now that the original technical difficulties have been overcome, 
we invert this chronological order. We start from the most general case, i.e., the 
massive GN model at finite temperature and chemical potential, and specialize further 
wherever additional analytical results or physics insights are available. In this way we 
can hopefully convey a more coherent picture of what has been learned in the meantime. 

This paper is organized as follows. In section 2, we remind the reader of the 
state of the art around the year 2000 and recall our criticism on the widely accepted 
phase diagram of the GN model at that time. We motivate the analytical form of 
the correct self-consistent HF potential in section 3 and sketch the calculation of the 
grand canonical potential and the proof of self-consistency in section 4. Section 5 
exhibits the revised phase diagram "in full glory". In the following four sections we 
have added various results in the form of figures and analytical formulas which have 
been obtained in simpler special cases, letting one or both of the parameters T and 7 go 
to zero. If we take the low density limit p — >^ at zero temperature, we also recover the 
properties of individual baryons as summarized in section 10. Section 11 contains a brief 
overview of the many cross-relations between the GN model as a relativistic QFT and 
the theory of certain quasi-one-dimensional condensed matter systems. It is followed by 
our concluding remarks in section 12. 

2. Reminder of the situation five years ago 

Let us briefiy recall the phase diagram of the GN model as it had been widely accepted 
by the year 2000. This is necessary to put the recent development into perspective. 
Besides, we will see that the results remain valid in some regions of the phase diagram, 
including some critical curves. The most basic quantity at finite temperature T and 
chemical potential fi is the grand canonical potential. In the large N limit, it can be 
computed either canonically via the relativistic HF approach or by using semi-classical 
functional integral methods. In the more popular path integral approach, one introduces 
an auxiliary bosonic field a and integrates out the fermions exactly (Gaussian integral 
over Grassmann variables), then applies a saddle point approximation to the remaining 
bosonic functional integral. In the original analysis at 7 = fT] and its generalization 
to finite bare fermion masses ^21; it has been assumed that the classical value of a (i.e., 
the saddle point) is spatially constant, yielding a (/i, T)-dependent dynamical fermion 
mass M. This mass is found by minimizing the renormalized grand canonical potential 
density 

* ^ ^ D - (1^ - v) r ^^'4(^ — )] w 

{l3 = l/TjUi = 1, E = a/^^ + M'^ ) with respect to M. Equivalently, one can solve the 
self-consistency condition for the (thermal) fermion condensate, 

M = mo-Ng^{^ij),^. (5) 
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Depending on the parameters (7,/i,T), \1/ possesses one or two local minima with the 
possibility of a first order phase transition. The phase diagram in the chiral limit 7 = is 
shown in figure 1, whereas the 3d-plot in figure 2 exhibits the phase structure in the full 
(7, /X, T)-space. Let us first look at figure 1. The line AB is a critical line of second order 
transitions (the thermodynamic potential changes from one minimum to a maximum 
and a minimum). The point B is a tricritical point located at = 0.3183, fi = 0.6082 
(all numbers here are in units of m); it separates the second order line from a first 
order line BD along which the potential has two degenerate minima and a maximum. 
The endpoint D lies at /i = 1/V2 where the T = phase transition occurs. Lines BC 
and BE are boundaries of metastability; when crossed, the potential acquires or loses 
a second minimum. In region OABD, chiral symmetry is broken and the fermions are 
massive; the outside region has unbroken chiral symmetry and massless fermions. As 
the parameter 7 is switched on, the 2nd order line AB disappears in favor of a cross-over 
where the fermion mass changes rapidly, but smoothly. The first order line on the other 
hand survives, ending at a critical point. If plotted against 7, these critical points lie 
on the third curve PtC emanating from the tricritical point as shown in figure 2. For 
7 > the effective fermion mass M is everywhere different from zero. If one crosses 
the shaded critical "sheet" in figure 2, the mass changes discontinuously, dropping with 
increasing chemical potential. 
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Figure 1. "Old" phase diagram of the massless GN model, assuming unbroken 
translational symmetry. For a discussion, see main text. Units of m; adapted from 
reference 



These results looked rather convincing and indeed have been confirmed repeatedly 
by other authors, see e.g. [13 CHI- Nevertheless, they suffer from one disease which was 
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Figure 2. "Old" phase diagram of the massive GN model as a function of 7,^,T, 
assuming unbroken translational symmetry ^| I14L 1^ . The various phase boundaries 
are explained in the main text. The 7 = plane corresponds to figure 1, the tricritical 
point Pt to the point B there. 



pointed out in reference The GN model possesses baryons (multi-fermion bound 
states) due to binding effects which are not suppressed. Consider the low density 
limit of baryonic matter at T = 0. Since the baryon-baryon interaction is known to be 
short-ranged and repulsive, one expects widely spaced baryons in the low density limit, 
implying the following slope of the energy density at the origin, 

d£ 



dp 



= Mb. (6) 

p=0 

Here, Mb is the baryon mass. Conversely, since equation © is just the chemical 
potential at T = 0, the phase transition at T = should occur at the critical chemical 
potential fi = Mb- The phase diagrams in figures 1 and 2 are in conflict with this 
expectation. At 7 = for instance, the first order phase transition takes place at 
/i = l/\/2, whereas the baryon mass is Mb = 2/11 (always taking out a trivial factor of 
A^). At the time of writing the article it was already clear that the problem had to 
do with the assumption of translational invariance. Indeed, if the single baryon breaks 
translational invariance, there is no good reason why this symmetry should not be broken 
at finite baryon densities as well. (One can actually understand the value of of 
the critical chemical potential in terms of a droplet model for baryons, characteristic for 
the mixed phase at low density. The argument which was first developed in the context 
of the GN model with continuous chiral symmetry applies here as well.) 
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Incidentally, a related observation had already been made in 1987 by Karsch et 
al. whose numerical lattice calculation gave first hints that kink-antikink configurations 
play a role in the phase diagram at low density fHI- However the authors concluded 
erroneously that there was a problem with the mean field approach, an opinion which 
has often been voiced in the literature since then. This calculation was too crude to 
determine the phase diagram quantitatively, in particular the statements about the 
order of the phase transitions are no longer tenable, as we now know. 

Unfortunately, it is not easy to relax the assumption of a constant potential, 
since this requires a major new effort. This can already be judged from the rather 
involved single baryon problem |H] which must be contained as a limit in the full 
calculation. Which potential S{x) should replace the constant mass M so as to cure the 
aforementioned discrepancy? This will be answered in the next section. 

3. How to choose the right ansatz for the scalar potential 

Our strategy within the relativistic HF approach is extremely simple in principle: guess 
the functional form of the HF potential and verify its self-consistency. At finite density 
we expect the appearance of a regular array of baryons, hence we will be searching for 
periodic solutions. But what guidance do we have to pick the right ansatz? Is there any 
reason at all to assume that the exact potential can be given in closed analytical form? 
The following heuristic considerations should be of some help. 

Owing to the (-0-0) ^-interaction in the Lagrangian (^, the Dirac-HF equation for 
the GN model assumes the form 



of the 7-matrices has the advantage that the equations for the upper and lower 
components (j)± of the Dirac spinor ip can be decoupled by simply applying the Dirac 
Hamiltonian twice. 



Equation Q states that the Schrodinger-type Hamiltonians with potentials U± = S'^±S' 
have the same spectra, a textbook example of supersymmetric quantum mechanics with 
superpotential S. S{x) in turn depends on the eigenfunctions ipa and eigenvalues Ea 
through a self-consistency relation. 




(7) 



with real scalar potential S{x). The representation 

7° = -cTi, 7^ = io-3, 7^ = 7°7^ = -(T2 



(8) 




(9) 




1 



1 



(10) 



It generahzes equation (jSJ to an x-dependent scalar potential. 
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S^±S' = - ^ , , ^ (12) 



Let us first recall the results for single baryons in the massive GN model |19t I2(jj. 
The scalar potential for a baryon has the form 

S{x) = I + y [tanh.{yx — cq) — tanh(?/x + cq)] , ~ 2 ^^^^'^^^^ ^^^^ 

The parameter y depends on the bare fermion mass and the number of valence fermions. 
Here the corresponding isospectral potentials U± in the 2nd order equation Q are given 
by the simplest Poschl- Teller potential 

cosh^(yx ± Co) 

and differ only by a translation in space. The distinguishing feature of potential ()12|) is 
the fact that it is refiectionless; indeed, this is the unique reflect ionless potential with 
a single bound state. It is well known that static solutions of the GN model must 
correspond to refiectionless Schrodinger potentials |HII221- The fact that the ansatz (|TT|l 
leads to self-consistency is therefore quite plausible. 

Take now a lattice of infinitely many, equidistant Poschl- Teller potential wells. As 
pointed out in references [23123) ^^e lattice sum can be performed yielding a Lame-type 
potential, 

E ; = f— V I ^ - (-■^) I ■ (13) 

„iHoCosh2(a;-nd) V ^ / \^^' V ^ // 

(K, E denote complete elliptic integrals of first and second kind, K', E' the 
complementary ones with argument k' = \^1 — k?, sn a Jacobi elliptic function of 
modulus k) . Comparing the spatial period of both sides of equation (fT!^ , we can relate 
d and k via 

K , , 

d = n—. (14) 



How does the fact that the single potential wells ()12|1 are refiectionless manifest itself 
in the periodic extension (jl3j) ? This has been answered long time ago [2^1 12^]: The 
periodic potential has a single gap (or, in general, a finite number of gaps), in contrast 
to generic periodic potentials with infinitely many gaps. Thus refiectionless potentials 
generalize to "finite band potentials" as one proceeds from a single well to a periodic 
array. In the same way as the sech^-potential is the unique refiectionless potential with 
one bound state, the sn^-potential is the unique single band potential. Guided by these 
considerations, let us try to find the most general superpotential of the (single gap) 
Lame potential. After a scale transformation 

S{x) = AS{0 ^ = Ax, (15) 

U±{x) = A'^U±{^) should assume the form of the Lame potential plus constant, up to a 
possible shift of the rescaled coordinate (^+ = ^ + b), 

= 52 + 5' = 2KWe+ + r/, 
U_ = - S' = 2kV^ + 7], (16) 
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or, equivalently, 

= (sn^e^ + sn^e) + V, 
S' = (sn^e+ - sn^e) . (17) 

Let us try to solve equations (fTTj) for 5* and rj. We differentiate the upper equation (fTTjl 
using 

(snO' = cn^dn^ (18) 

and divide the result by the lower equation, obtaining 
~ sn^+cn^+dn^+ + sn^cn^dn^ 

sn^,^+ — sn^4 

By specializing equation |T7jl to ^ = we can also determine the constant r], 



r,= \-l-t,\ (20) 
sn^o 

Incidentally, S can be cast into the somewhat simpler form 

;=,/^N 9 , ^ ,\ cn6dn6 

5(0 = sn6sn^sn(^ + 6) + — . (21) 

sn 

The scale factor A in equation ()15|) is not constrained by these considerations, so that 
the final answer for S{x) depends on three real parameters A, n, b. Thus we conclude 
that the potential P5I19|) is the most general Dirac potential leading to a single gap 
Lame potential (plus constant) in the corresponding 2nd order equations. This explains 
why it is a good starting point for finding periodic, static solutions, and has turned out 
to be the key to the phase diagram of the massive GN model. Incidentally, we did not 
derive the ansatz for S{x) in this way originally, but took it over from a mathematically 
closely related problem which had already been solved in condensed matter physics, the 
bipolaron crystal in non-degenerate conducting polymers. We will have more to say 
about this relationship in section 11. Notice also that the simple and intuitive relation 
between the single baryon and the crystal exhibited in equation (fT^ is deeply hidden 
in the corresponding Dirac potentials pijl and p5|19j) due to the non-linear, non-local 
relationship between S and U±. 

Inserting our ansatz into equation (jH)), we arrive by construction at the single gap 
Lame equation in the form 

(-^ + 2«:'sn2(e +{bT b)/2)j 0± = S(jy±. (22) 

Using equation ^lU^. the relation between Dirac eigenvalues E and Lame eigenvalues £ 
is 

^ = 5 - ^ + ^ + (23) 
A'^ sn'^b 

The solutions of equation ()22j] are well known since reference |2ZI. The eigenfunctions 
are fairly complicated, although they can still be given in closed form, 

<^+(0=A/'^^|(|f^e-(")^. (24) 
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Here, H, G and Z are the Jacobi eta, theta and zeta function, respectively. The 
parameter a is related to the energy and Bloch momentum. An example of the dispersion 
relation is shown in figure 3. 




upper band 



lower band 



Figure 3. The dispersion relation of the Lame potential for k — 0.8, showing a single 
gap; to = E/A and p = q/A are (reduced) energy and Bloch momentum, respectively 
(from reference |28|V 



The fact that S{x) is a finite band potential is important for getting self-consistency. 
At the same time, it makes the problem analytically tractable. Of course, at this point 
one cannot be sure that the simplest finite band potential is sufficient to solve the massive 
GN model at all temperatures and chemical potential, but this is what eventually will 
come out. 



4. Minimizing the grand potential and self-consistency condition 

Having picked the ansatz for the scalar potential S{x), we now have to minimize the 
grand potential with respect to the three parameters A, k, b and verify that it implies self- 
consistency of the ^/'^/'-condensate. In HF approximation, the grand canonical potential 
density per flavor consists of the independent particle contribution and the double 
counting correction to the interaction \l/2. 



^1 
^5 



^1 + ^2, 

1 r^/'^ 



— / dg In (l + 



,/3(£-f/.) 



2Ng^ 



dx{S{x) — mo) 



1 + e' 



2K/A 



(25) 
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(in the path integral approach, and \E'2 correspond to the fermion determinant and 

the tree level term, respectively). A/2 is an UV cutoff which will eventually be sent to 

infinity, i the spatial period of S{x), E the single particle energy. 

In \l/2, we need the spatial averages of 5* and S"^. It is convenient to introduce three 

basic functions of b and k, 

1 cn6dn6 E 
s=^-, t = -— . u = l-—, (26) 

where s and t are related by 

t^ = s{s-l){s-K^). (27) 
The spatial averages can then be written in the compact form 

{S) =AiZ + t/s), 

(5') = (s - 1 + 2m - k") (28) 

with Jacobi's Zeta function Z = Z(6, k). 

Turning to ^i, let us first transform the momentum integral into an integral over 
Dirac energies. The density of states for the Lame potential [221 implies a change of 
integration measure {p = q/A,uj = E/A) 

dp uj{ij'^ — s + u) 

W ={uj^-s + l)(cu2 - s + K^){uj^ - s), (29) 

where the plus sign refers to the upper band, the minus sign to the lower band. With 
the shorthand notation 

a = (3A, v = fi(3, (30) 

can be written as the following integral over the allowed bands 

7r/?2^i = -a (^j^^duj + j^^<^^ ^ In [(1 + e-""+'^)(l + e'^^")] . (31) 

The energy cutoff A^ has to be computed from the momentum cutoff A/2 and equation 

A. = ^ + ^ + 0(A-^). (32) 

Due to the quadratic divergence of the integral, it is necessary to keep the next-to- 
leading order term here. Following references |28t l30j . we combine the integral over 
both energy bands as well as over positive and negative energy modes into the real part 
of a line integral in the complex cj-plane. The path of integration runs infinitesimally 
above or below the real axis, 

7r/?2^i = -a lim Re ^^ ^(^'-^+») ^ (33^ 



We have to minimize \1/ with respect to A, 6, k for fixed /x, /3, 7, a rather tedious task at 
first glance. In order to simplify the computations, we use the freedom to minimize \& 
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with respect to any other set of independent variables. Specifically, we propose to replace 
A, b, K by the spatial averages (S*), (S*^) and the spatial period £ of S, see equations 
and (j2Hl)- It turns out that these variables are extremely convenient for the proof of 
self-consistency. The stationarity conditions for the grand potential then read 



d{S) 

d 

d{S^) 



PFo = 
P^Fi/2 



0, 



0, 



0, 



(34) 



where we have introduced functions Fi which play a key role in our approach. The 
derivatives with respect to the new, composite variables can be taken trivially in the 
case of \l/2 since 

1 



2Ng'' 



[{S^)-2mo{S)+ml 



(35) 



Unfortunately this is not true for available only in terms of the original variables 
A, b, K in equation (jH^ . In order to compute F^, we therefore invoke the chain rule, 



dA 
9. 



(36) 



dA{S) dA{S') dAi 

[ d,{S) d^{S') dj J 

Upon evaluating and inverting the Jacobian matrix on the right hand side of equation 
(IHUj). the Fi can be expressed in terms of directly computable derivatives of After 
some algebra [301 we arrive at 



PFo 
P'F,/2 
V a^F2/£ J 



Fo 



Fi 




(37) 



All three functions Fi vanish at the minimum of the thermodynamic potential ()25p. 

Equations (jSBI), and (jHTjl are the basis for our computations. They are not 
yet in a form suitable for numerical calculations since they still involve bare parameters 
and a cutoff, but the present form is illuminating with respect to self-consistency. We 
therefore first turn to the issue of self-consistency. 

The self-consistency condition for the scalar potential S{x) was given in equation 
(Hn)). For a single mode, '^'^ reads [HI] 
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We transform the thermal expectation value of ipip once again into a complex integral, 

1 /■A/2 - / 1 1 \ 

A , f^+i^ , uj'^S-t 1 , , 

= -Re lim / du — ^ . 39 



A partial integration (picking up a boundary term) then enables us to express the 
thermal expectation value in terms of the functions Fo, Fi introduced above, 

In the minimum of the grand canonical potential, Fo and Fi vanish. This reduces 
equation pn|) to equation dHH), thereby establishing exact self-consistency of the scalar 
condensate. 

We now comment on the practical procedure to determine the phase diagram and 
thermodynamics of the massive GN model within this framework. First, we have to 
eliminate the bare parameters Ng"^, mo and the cutoff A in the standard way. Since all 
UV divergences are due to vacuum effects, there are no new difficulties as compared to 
the T = case. All we need is the vacuum gap equation 

^ ^ ;i + mo)lnA= -(lnA + 7). (41) 



Ng"^ TT 71 
The dynamical fermion mass in the vacuum (set equal to 1) and the confinement 
parameter 7 are physical parameters which are kept fixed while one lets Ng'^ —>■ 0, 
mo — *■ and A — >■ cxo. Irrelevant divergent terms — A^/Stt and — /iA/27r from the energy 
and baryon density of the Dirac sea can simply be dropped. An expression for the 
renormalized grand canonical potential in which the limit A ^ 00 can safely be taken 
is 



7r/?'^ren= lim 

A^oo 



+ \ (In A + 7 - 1) - ap-f{S) 



— aRe lim / du — In (2 cosh 



^ , . (42) 

In principle, we have to solve the three equations Fi = simultaneously. It is 
worthwhile to examine more closely the way in which the Fi depend on the 6 relevant 
parameters (a, b, k, 7, u, (3), since this suggests a simpler strategy of how to minimize 
We first note that F2 is a convergent integral which does not require any regularization 
or renormalization. Since F2 depends neither on 7 nor on (3, the equation 

F2{a,b,K,u) =0 (43) 

can be solved for a, say, for given b, k, v. Now let us focus on the 7, (3 dependence of the 
other two equations. Fo is also free of divergences. All we have to do here is to replace 
the ratio of bare parameters Timo/Ng"^ by 7. The equation Fq = is then turned into 

•jP = atIo{a,b, K,!^), (44) 
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where Jq can be inferred from equation (jSTj). The integral in Fi on the other hand has 
a logarithmic divergence at the lower integration limit. Isolating the divergence and 
eliminating l/Ng"^ with the help of the gap equation ()4H1 . the In A terms are cancelled 
and the equation Fi = assumes the form 

7 - In/? = /i(a, 6, z/). (45) 

Equations and (jlKjl can be combined into 

7 + In 7 = Ji + In (at/o) , 

(3 = ath/i. (46) 

The first of these equations can be solved for 7, the second one then yields an explicit 
expression for (3. In total we have reduced the problem of finding the minimum of a 
function of three variables to the simpler problem of finding the zero's of two functions 
of one variable each. 

Throughout this section, we have made use of complex integration so as to 
exhibit the formal structure of various expressions and the self-consistency in the most 
transparent way. In order to actually compute the one-dimensional numerical integrals 
in ^ and in the Fj, it is advisable to convert the integrals back to standard real integrals 
over the allowed energy bands. For more technical details and explicit expressions for 
Jo, /i, F2 we refer the reader to reference 



5. The full phase diagram of the Gross-Neveu model 

In the previous section, we have outlined how to minimize the grand canonical potential 
given the 3-parameter ansatz for S{x). This can be used to evaluate bulk thermodynamic 
observables for any given T, /i and 7. In the present section we turn to the structure 
of the phase diagram, including the location of the phase boundaries, the order of the 
transitions and the symmetries of the various phases. The following considerations have 
proven useful for this purpose. We expect phase transitions between a crystal phase and 
a homogeneous phase where S{x) = M is constant. In the chiral limit, we can further 
distinguish between a chirally restored phase (M = 0) and a massive, chirally broken 
phase (M > 0). We therefore first need to understand for which choice of parameters 
S'(x), equations ()15I19|) . goes over into a constant. The elliptic modulus n varies between 
and 1. At the two boundaries of this interval, we find 

lim 5* = cot 6, 

lim S = cothfe + tanh^ - tanh(^ + b). (47) 

The first limit is a constant, the second describes the single baryon profile. Clearly, both 
limits are relevant for phase boundaries where S goes over from a periodic function in 
the crystal phase to a constant one in the massive Fermi gas phase (M = A cot b or 
Acothb, respectively). As far as b is concerned, important special values are 6 = and 
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The value 6 = is the only value of b where k = and k = 1 can coexist and therefore 
plays a prominent role in the phase diagram. In the chiral limit, b takes on the value K. 
If we let 6 — > K and k — simultaneously, S vanishes and we can connect the crystal 
phase continuously to a chirally restored phase. 

These observations are the key for computing the phase boundaries separating 
the crystal phase from the Fermi gas. Evidently, they are only applicable if the phase 
transitions are continuous, but this is indeed what we find. For k = or 1, the functions 
Fi further simplify, so that the actual computation along the lines described in section 4 
is quite manageable. We have computed lines of constant b and constant u, since this is 
most easily done in our scheme. For each value of k, a certain (6, z/)-grid is mapped onto 
a two-dimensional surface in (7, /i, T)-space by minimizing the grand potential. The 
resulting curved surface represents a 2nd order phase boundary. 




Figure 4. Phase diagram of the GN model in the chiral limit |28[ ■ All phase 
boundaries correspond to second order transitions. 



We now turn to the numerical results, starting with the chiral limit 7 = and 
focusing on the phase boundaries. These are depicted in figure 4 which supersedes 
figure 1. The 2nd order line AB of the old phase diagram is unaffected. The first order 
line BD of figure 1 is replaced by two 2nd order lines delimiting a novel kink-antikink 
crystal phase. The tricritical point B is turned into another kind of multicritical point 
labeled Pl in figure 4, located at precisely the same (/i, T) values. As we switch on 7 
(figure 5), the second order line separating massive (M > 0) and massless [M = 0) 
phases disappears as a consequence of the explicit breaking of chiral symmetry. The 
crystal phase survives at all values of 7, but is confined to decreasing temperatures 
with increasing 7. For fixed 7, it is bounded by two 2nd order lines joining in a cusp. 
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The cusp coincides with the critical point of the old phase diagram but has once again 
a significantly different character. The crystal phase exists and is thermodynamically 
stable inside the tent-like structure formed out of two sheets denoted as / and //. These 
sheets are defined by k = (/) and k = 1 {II), respectively. The line PtC where they 
join corresponds to 6 = and coincides with line PtC in figure 2. The baseline of sheet // 
in the [fi, 7)-plane has a simple physical interpretation: It refiects the 7-dependence of 
the baryon mass in the massive GN model, or equivalently the critical chemical potential 
at T = 0. The chiral limit 7 — > can be identified with 6 — K. 

Another view of the phase diagram is displayed in figure 6. Here we plot the phase 
boundaries in the {fi, T)-plane for several values of 7. These two-dimensional graphs 
correspond to cutting the three-dimensional graph in figure 5 by (equidistant) planes 
7 = constant. They provide a better view of how the two critical lines are joined in 
a cusp and exhibit that the region in which the crystal is thermodynamically stable 
shrinks with increasing bare fermion mass. 

The cusp line where the sheets / and // are glued together can be determined 
by taking the limit b ^ in Fi. As discussed above, in this limit the scalar 
potential becomes both homogeneous and k- independent. One can show analytically 
that the resulting equations Fj = for i = 0,1,2 are equivalent to the conditions 
vl>' = 0, \E'" = 0,\E^"' = (' = Om) for the translationally invariant calculation 
Hence the curve 6 = in the new phase diagram coincides with the line of critical 
points of the old phase diagram. One can also convince oneself that the critical first 
order phase transition sheet shown in figure 2 is tangential to both sheets / and // along 
the line of endpoints. 

6. Ginzburg-Landau theory near the tricritical point 

In the vicinity of the tricritical point Ft, we can derive a Ginzburg-Landau effective 
action by expanding ^'ren in a = pA |2Hl IHOj- Since S{x) is both weak and slowly 
varying there, this could also be obtained approximately without knowing the exact 
solution. Here, we turn things around and derive it via a Taylor expansion from the full 
expression. This should help to understand the character of the multicritical point. 
We therefore start with an expansion of \E'ren in equation (jl^ in powers of a. 



Here, ip{n, z) is the polygamma function, the ra-th derivative of the digamma function 
ip{z) = r' {z)/r{z). The coefficients Cn can be obtained from a generating function. 




(49) 
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Figure 5. Revised phase diagram of the massive Gross-Neveu model jHlEOl- The 
shaded surfaces /, // separate the kink-antikink crystal from a massive Fermi gas and 
correspond to 2nd order phase transitions. The 7 = plane is the same as figure 4. 



One can now verify that the combinations c^+i + (m — s)c„ appearing in equation 
are related as follows to spatial averages of powers of S and its derivatives with respect 
to X, 

1 



ci + {u- s)co 

C2 + {u- s)ci 
C3 + (m - S)C2 



2^2 
3 

8A^ 
5 



16^6 



{{s') + {{ST)] 

{{s')+l{{sr)+5{s\sr)). 



This enables us to write down the Ginzburg-Landau effective action in the form 

^efr(7) = ^eff(7 = 0) + ^ - 2S) , 
where the effective action at 7 = is 



(51) 
(52) 
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56 + \s"f + bS\S'f\ Re ^ f4, - + \] 



2-kTJ 



(53) 



The instability with respect to crystallization is related to the fact that the "kinetic" 
term ~ [S'Y can change sign, depending on ^ and T. We will come back to this formula 
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Figure 6. Two-dimensional sections 7 — const, through the phase diagram of figure 5. 
The fat hues belong to 7 = 0, the thin lines to 7 = 0.1, 0.2, 1.2, from top to bottom. 
The position of the cusp agrees with the critical point in the old phase diagram. From 
reference PH] . 

in section 11 where we point out that it has a hteral correspondence in condensed matter 
physics, namely in the theory of inhomogeneous superconductors. 

Notice that if one drops all derivatives of S in the effective action, one gets the 
result for the GN model under the assumption of unbroken translational invariance. 
The tricritical point for instance is defined by the simultaneous vanishing of the 5*^ and 
coefficients which happens at fit = 0.608221, Tt = 0.318329. Since the coefficients 
are the same as in the full effective action, we can understand why the tricritical point 
stays at the same place even though we allow for x-dependent potentials. 



7. More about the chiral limit 



In references |2H1IS21, a number of additional results has been obtained in the chiral limit 
7 = 0. Here we give a selection of figures and formulas to highlight certain features of 
this simpler special case. 

We first reiterate that only at 7 = we are dealing with three distinct phases: the 
crystal, a massive and a massless Fermi gas. The scalar potential can be obtained by 
letting 6 — >• K in equations ()15fl9j) . 

SiO = = '^'snesn(e + K) (54) 
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and has a higher symmetry than for 7 7^ 0, namely 

S{x + e/2) = -S{x) (55) 

where £ is the spatial period. This is actually a remnant of the original discrete chiral 
symmetry of the model. Translational invariance and the 7^ transformation both break 
down, leaving the unbroken discrete symmetry ()55|) . i.e., a translation by half a period 
combined with a 7^ transformation. 

An expanded plot which reveals more details about the shape of the phase boundary 
separating crystal and massive Fermi gas (labeled k = 1) is displayed in figure 7. 
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Figure 7. Close-up on the phase boundary between homogeneous and inhomogeneous 
ordered phases (k = 1). Notice the different scale on the /[^-axis as compared to figure 
4. From reference |28| . 

A way of presenting the phase diagram complementary to that in figure 4 is given in 
figure 8 where we have transformed all phase boundaries from the (T, /i)- into the (T, p)- 
plane. Here, the first order line of the old solution splits up into the two dashed lines 
which delimit the mixed phase region (droplets of chirally restored matter in the chirally 
broken vacuum). This should be replaced now by the two solid lines going downward 
from the tricritical point and enclosing the crystal phase. At T = in particular, the 
crystal phase is stable at all densities. 

As we approach the k = phase boundary, S ~ cos2gx (with an amplitude 
vanishing at k = 0), and the wave number q can serve as order parameter for the 
breakdown of translational invariance. In figure 9 we show the dependence of this order 
parameter on fi as one moves along the phase boundary. The solid line is the curve q = p 
which is approached asymptotically by the full calculation. At /i = fit, the tricritical 
point of the old solution, we see a clear signal of a 2nd order phase transition with 
breakdown of translational invariance. 
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Figure 8. Revised phase diagram of the GN model in the (T, p)-plane (32] ■ The 
dashed hnes belong to the old phase diagram, where they enclose the mixed phase. 
This "droplet" region is superseded by the crystal phase featuring baryons. 




fJ-t 



Figure 9. Wave number characterizing the crystal period along the phase boundary, 
showing a continuous phase transition at the tricritical point /i^. From reference 



Since the phase boundaries are of particular interest and easier to compute than 
other thermodynamic quantities, let us give the explicit formulas derived in reference 
[2Hj from which they can be obtained. 

The boundary between the chirally restored phase and the crystal is characterized 
by K = 0. It can either be derived by almost degenerate perturbation theory |32j or 
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from the full thermodynamic potential. In the latter straightforward expansion 

around k = yields 



In = - min Re 

AttT 2 a>0 



(56) 



The resulting curve is labeled "k = 0" in figures 4 and 7. The asymptotic behavior for 
large /x can also be determined, 

^- = 5^- 

The phase boundary does not reach the T = axis so that chiral symmetry is not 
restored at T = no matter how high the density is, at variance with the naive 
expectation for an asymptotically free theory. 

For small u where Re ip (2,^ + < 0, the unique minimum is at a = 0. In this 
range of u the phase boundary does not touch the crystal region in the (/i, T)-diagram. 
It corresponds to the transition between the massless and the massive homogeneous 
solutions described by 

Upon using ip{l/2) = —C — ln4 (C ~ 0.5772 is the Euler constant), we reproduce the 
well-known value of the critical temperature at /i = 0, 

e^ 

= — . (59) 

TT 

Next we turn to the non-perturbative phase boundary at k ^ 1. It separates the 
crystal from the massive phase. The relation between a and u along the phase boundary 
can be deduced from 



d 1 ""r^ 1 r (i + ^(z/ + acosv?), 

/ d(^^Imln— >^ ^-^ ^ = 0. (60) 

daa J cosip r (^^ + ^{u - acosip) 



With these values for a and z/, 



, /5 ^ 1 
- In h - 



/ dipRe^(l + ^{u + acosip)) =0 (61) 
Jo \2 2tt / 

yields /3 along the phase boundary. The resulting curve is shown in the phase diagram 

in figures 4 and 7 (label "k = 1"). Equation ()6H) alone with values of (a,!/) that are 

not restricted by equation (jHIHl gives the connection between a (or, equivalently, the 

effective fermion mass), u and (3 in the massive phase outside the crystal region. 



8. More about the zero temperature limit 

Let us begin with some additional plots on cold, dense matter in the massive GN model 
p?T] . First of all we wish to illustrate how the self-consistent scalar potential depends on 
7. This is exhibited in figure 10 at low and high density, respectively. The deepest curves 
always correspond to 7 = 0, where the potential oscillates symmetrically around zero 
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Figure 10. Self-consistent scalar potential S{x) versus x for pf = 0.3 (left) and 
Pf = 2.5 (right). From bottom to top: 7 — 0,0.01,0.1,0.3,0.75,2.3. From reference 

m 



due to the residual symmetry As we increase the symmetry violation parameter 

7, the potential oscillates with decreasing amplitude around a value close to the mass 
M which the fermions would acquire in the translationally invariant solution. It is 
surprising that such a variety of potential shapes in the Dirac equation can all be 
reduced to the standard single gap Lame equation. 

The next result which we should like to show is how the density varies with the 
chemical potential. At T = 0, the chemical potential can be obtained by differentiating 
the energy density with respect to the mean fermion density, 

Pf 

op n 

If we assume unbroken translational invariance (figure 11, left), we find discontinuities in 
these curves, confirming the result of reference jT2| about a first order phase transition. 
Repeating the same calculation for the crystal solution (which is the stable one), all 
the curves become continuous, signaling a 2nd order phase transition (figure 11, right). 
The critical chemical potential in this latter case coincides with the baryon mass, as 
expected on general grounds. By contrast, the first order transition in figure 6 happens 
at a chemical potential which has at best the meaning of an approximate baryon mass 
in a kind of droplet model, cf. the discussion in section 2. 

The two real parameters which determine S{x) at T = are the elliptic modulus 
K and the shift b. The scale factor A on the other hand is determined by the baryon 
density and k, 

A = MS. (63) 

71 
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Figure 11. Baryon density versus chemical potential [HI]- Left: (unstable) 
translationally invariant solution, showing a first order transition. Right: (stable) 
crystal solution with 2nd order transition. Curves ft'om left to right: 7 = 
0.01,0.1,0.3,0.75,2.3. 



Expression implies that the mean density scales with the inverse spatial period of 
the potential, a consequence of the fact that the valence band is completely filled (for 
matter) or empty (for antimatter). 

We now turn to additional analytic results. At T = 0, all integrals can be reduced 
to incomplete elliptic integrals of first, second and third type F, E, U. We illustrate this 
fact in the case of the ground state energy and the self-consistency conditions. The 
ground state energy density is split up according to Eg ^. = Ei + E2 where Ei is the sum 
over single particle energies of the occupied states, E2 the double counting correction 
to the interaction energy. For Ei we get 
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+ xE{p, q) 



(64) 



- n (p, 

with the following definitions [t] has been defined in equation (j20|l ]. 




F(p, q). 



p = dn6, 

q = k' /dnb. (65) 

E2 on the other hand coincides with \l/2 in equation ()25j) since the HF double counting 
correction is T-independent. After renormalization and upon using the vacuum gap 
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equation together with equations (PH|) . one finds 

A2 A 

^2 = — (7 + In A) (s - 1 + 2u - K^) - -7(Z + t/s). (66) 

When adding up Ei and i?2 the logarithmic divergence drops out and a finite result 
depending only on physical parameters is obtained. 

Consider the self-consistency condition next. For 7 7^ 0, it can be cast into the 

form 

= Acnb F(p, q) — 'J sn^b, (67) 



J = Acnb 



n (p, {k'Y, q)+x{l + HA^/2 + ^)^ 



The original x-dependent condition has been transformed into two x-independent 
equations. They determine k and b for given pf and 7. 

We finish this section with a few words about the high density limit where the 
calculation becomes perturbative (see jHI])- We use almost degenerate perturbation 
theory, following the standard weak coupling approach from solid state physics 
textbooks. In the double-counting correction, we have to take into account the bare 
fermion mass, and we allow for S'o 7^ in addition to S±i 7^ [Si are the Fourier 
components of the periodic potential S{x)]. Thus our present ansatz for S{x) is 

S{x) = So + 2Si cos{2pfx). (68) 

The approximate energy density becomes 

Minimizing with respect to Sq and 5*1, we find 
5'o[ln(2p;)+7]-7 = 0, 

^i[27 + ln(4p^5i)] = 0. (70) 
The first equation has the unique solution 

^° = 7 + ln(2p,)- ^''^ 
The 2nd equation has two solutions: 5*1 = 0, corresponding to unbroken translational 
invariance, and 

S, = -^e-2^ (72) 
4pf 

for the soliton crystal. Comparing the energy densities of these two solutions, 

E,...{S, ^ 0) - E,,,XS, = 0) = -^e-^^ (73) 

we learn that the crystal is favored, but the energy difference decreases rapidly with 
increasing 7. We have shown this simple calculation in detail since it proves non- 
restoration of translational invariance at high density for arbitrary bare fermion mass. 
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9. More about the double limit T = 0,7 = 

The simplest problem beyond individual baryons and the first one which could be solved 
analytically |33j is the ground state of baryonic matter in the chiral limit. The ansatz 
for the reduced HF potential has the form (j54|l . whereas the scale factor A is given 
by equation (jHS))- Thus only a single variational parameter k is left. As expected, all 
analytical formulas of the previous section simplify tremendously in the chiral limit. 
Thus for instance, the renormalized ground state energy density can now be expressed 
entirely in terms of complete elliptic integrals 

E..- = ^ (4E + (K^ - 2)K) + MIS (2E + - 2)K) In (^) . (74) 
Minimizing with respect to k, yields the transcendental equation 

' = 1 = 2^^ ^''^ 

and self-consistency can be ascertained rather easily in this case. Eliminating pj from 
Eg,s. with the help of relation (fTSjl . we arrive at the following parametric representation 
of the ground state energy as a function of density in terms of the parameter k, 



An TTK^ V K 
Pf _ 1 



(77) 



TT 2kK 

In the low- or high-density limits, it becomes possible to systematically resolve the 
transcendental equation (|77jl . 

^ ^ l-8e-/^/ + ^^i^^±^e-^-/^/. 



►0 Pf 
1 1 



K ^ ^ + ^. (78) 

For the energy as a function of density, one then finds 

E,.s. ^ ^-t^ + tJ-F- (79) 



Pf- 



2n 2Hp) 2^Hpj 



In the low-density limit, the three terms correspond to the vacuum energy density, the 
contribution from the baryon mass (~ pMs with Mb = 2/7t) and a term reflecting the 
repulsive baryon-baryon interaction. At high densities, we can identify the free massless 
Fermi gas piece, the leading perturbative correction [32] and the next term coming from 
higher order effects, suggesting fast convergence. 
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We end this section with a few remarks about the baryon density. In the present 
case, the exact density is x-dependent and given by 

It has the following high and low-density limits: At low density (k —* 1) we recover the 
result for a single baryon 

p{x,K^l)^- r^ + - r2r^-— (81) 

4 cosh X 4 cosh [x/k + K) 

At high density (k — > 0) the total baryon density approaches a constant, 

Pla^-'^^O)^^, (82) 
unlike the scalar potential which keeps oscillating around with wave number 2pf. 

10. A note on baryons 

Single baryons are contained in the low temperature, low density limit of the general 
formalism, at least those with completely filled or empty valence state. They have also 
been studied as such, first in the chiral limit jS] and more recently in the massive GN 
model [ini 1201 • One finds that the scalar potential has the same shape as the niQ = 
baryon with partial filling of the valence level. 



with 



S{x) = l + y (tanh ^_ - tanh ^+) (83) 

^± = yx ± -artanh?/, (84) 

where m = 1 and y G [0, 1] is the only variational parameter. It depends on the 
occupation of the pair of valence levels and 7. Evaluating the ground state energy Mb 
and varying with respect to y, 

dy 

one obtains 



0, (85) 



^ = ^ + 2tan^ (86) 

2 TT 7r 

where we have introduced the angle 6 via y = sin 6 {0 < 6 < 11/2). Here u is the number 
of valence particles minus the number of holes in the negative energy valence state. The 
baryon mass at the minimum becomes 
Mb 2 

— — = — [sin 9 + J artanh(sin 9)] . (87) 

N TT 

The phase boundary in the T = plane of figure 5 agrees perfectly with the function 
Ms (7) evaluated from equations (j86l87|) for u = 1. This was not the case in the old phase 
diagram of figure 2, confirming once again that the revised results are now internally 
consistent. 
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11. Relation to condensed matter physics 



From the point of view of relativistic QFT, the GN model is viewed as a soluble model 
for strong interaction physics, exhibiting phenomena like asymptotic freedom, chiral 
symmetry breaking, dynamical mass generation, meson and baryon bound states. On 
the other hand, it can model quasi-one-dimensional condensed matter systems in the 
vicinity of a half-filled band. In this context, it has enjoyed success in explaining real 
experimental data, something unthinkable in particle physics due to its toy model 
character. The two most striking examples which we are aware of are conducting 
polymers on the one hand and (quasi-one-dimensional) inhomogeneous superconductors 
on the other hand. The basic physics is the same in both cases and closely related to the 
Peierls effect j33] of a one- dimensional electron-phonon system — dynamical formation 
of a gap at the Fermi surface resulting in a crystal structure. Details are different though, 
in particular the identification of the chemical potential is more subtle in the second 
case. Since the relationship between the GN model and condensed matter systems has 
played an important role for establishing the phase diagram of the GN model during 
the last few years, let us discuss it in somewhat more detail. 

1. The GN model and poly acetylene 

Conducting polymers have had a tremendous success story, culminating in the 2000 
Nobel Prize in chemistry for physicist A. J. Heeger and the chemists A. G. MacDiarmid 
and H. Shirakawa (for a review, see [SS])- A prominent example is trans-polyacetylene 
(PA). This polymer (CH)^; possesses two "dimerized", degenerate ground states with 
alternating short and long bonds (figure 12, left). Owing to a number of simplifying 
assumptions, its continuum description is mathematically equivalent to the symmetric 
(mo = 0) GN model, as was realized shortly after the seminal work of Su, Schrieffer and 
Heeger on the discrete model [SSI IHl IHZl IHH] • Dimerization plays the role of (discrete) 
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Figure 12. Left: Two degenerate, dimerized ground states a) and b) of frans-PA, 
described theoretically in terms of the massless GN model. Right: Two inequivalent 
configurations of czs-PA, leading to a description in terms of the massive GN model. 
Configuration c) has lower energy than d). 



chiral symmetry breaking in relativistic QFT. Solitons appear as kink-like defects where 



Updating the Gross-Neveu phase diagram 



28 



a transition between the two degenerate ground states happens, polarons are kink- 
antikink bound states. These excitations are important for understanding the electrical 
conductivity properties of doped PA (doping changes the number of electrons and is 
analogous to changing baryon number in the GN model). The solitons in particular 
have attracted considerable attention in the context of fermion number fractionization 
of which they are a prime example [^11 EHl 1^ • Solitons and polarons have originally 
been analyzed theoretically in close analogy with kink and kink-antikink baryons already 
known on the field theory side [H] . A good account of both the physics and the early 
history can be found in reference where also the limitations of the equivalence 
between the GN model and the theory of PA are critically examined. More recently, 
this theme has been taken up again by others Here, let us only mention a few facts 
which help to understand the curious equivalence of a non-relativistic electron-phonon 
system with a relativistic QFT. Particle/antiparticle degrees of freedom correspond to 
particles and holes, the half filled band to the Dirac sea. The role of the UV cutoff 
is taken over by the band width W (incidentally of the order of 10 eV in PA). The 
two spin-components (which exist only in the condensed matter case) are mapped onto 
2 flavors (i.e., the N = 2 GN model). In the widely used semi-classical approach, 
there is no difference between the condensed matter treatment and the large N limit 
used above. The linear ( "ultra-relativistic" ) dispersion relation for the massless fermions 
arises from a linearization at the Fermi surface, the Fermi velocity playing the role of the 
velocity of light. The lattice distortion (or phonon) field in the adiabatic approximation 
corresponds to the auxiliary scalar field a in the GN model, the gap parameter A{x) to 
the scalar potential S{x) and the electron-phonon coupling A to the coupling constant 
Ng"^. In the relativistic QFT case, one has to send the UV cutoff to infinity and the bare 
parameter Ng"^ to in a specific way, dictated by the renormalizability of the model. 
In polymer physics, W and A are physical observables to be taken from experiment. 
Nevertheless, in practice these differences do not matter for many questions. 

As mentioned above, the kink baryons are the solitons, the kink-antikink baryons 
the polarons, bipolarons and excitons of polymer physics. More specifically, for N = 2 
the valence level can only be empty, fully occupied or half occupied. In condensed 
matter physics, there is then a relation between charge and spin. For the questions 
discussed in the present paper, fully occupied or empty valence levels are most relevant 
— these should be identified with the bipolarons. In the massless GN model, the baryon 
with fully occupied valence level becomes a kink/antikink at infinite separation; this is 
the reason why bipolarons are not discussed in the case of degenerate polymers like 
trans-PA. 

If one starts to think about the massive GN model, it is not too hard to identify its 
condensed matter analogue: Conducting polymers with non-degenerate ground states 
like cis-PA (figure 12, right). Their theoretical description was initiated by Brazovskii 
and Kirova [44^ with the proposal that the gap parameter has two contributions, a 
constant, "external" one arising from the basic structure of the polymer and an x- 
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dependent, "internal" one due to electron-phonon coupling, 
A(x) = Ae + Ai(x). 



(88) 



If we identify A(a;) with the scalar mean field S{x) and Ae with the bare mass mo, we 
can immediately relate two problems from two different branches of physics. Baryons 
in the massive GN model then correspond to polarons and bipolarons in the polymer 

case nacni. 

Now consider periodic solutions of the gap parameter for both degenerate and 
non-degenerate polymers and compare them with relativistic QFT. In the degenerate 
polymer case, polaron crystal solutions to the continuum model (at T = 0) have been 
found and discussed in BSl They are mathematically identical to the ground 
state of baryonic matter in the chirally symmetric GN model, cf. section 9. In some 
cases the order parameter looks different from ours at first sight. However, the two 
expressions can be converted into each other by Landen's transformation for Jacobi 
functions in the form 



In the 80's and 90's, a lot of work was devoted to non-degenerate conducting polymers 
with the result that exact bipolaron lattice solutions were found by several groups. This 
suggests that the massive GN model should exhibit a soliton crystal at finite baryon 
density as well. It encouraged us to look for such solutions of the massive GN model 
and to try the functional form of the self-consistent scalar potential for non-degenerate 
polymers in the GN model. It has been derived by various methods such as inverse 
spectral theory |lHj, Poisson summation of periodic sums of single polarons jl9l 1501 l5T] 
or relation to the Toda lattice [23] — a rich source of inspiration also for those working on 
relativistic QFT. Reference [21] contains a particularly thorough discussion and corrects 
misprints in some of the original papers. In section 3 we have in fact presented the 
functional form determined in these works and used it as ansatz in a relativistic HF 
calculation. The arguments given in favor of this ansatz are also borrowed from the 
condensed matter literature. On the other hand we did not find the full phase diagram 
at finite T presented in section 5 in the polymer literature. 

2. The GN model and inhomogeneous superconductors 

Superconductivity is driven by fermion-fermion pairing (Cooper pairs), whereas 
the GN model features fermion-antifermion pairing (chiral symmetry breaking). We 
first have to understand how these two distinct physical phenomena are related to 
each other. Along the lines described in reference [221, one can actually map the GN 
Lagrangian onto a "dual" Lagrangian which has fermion-fermion pairing by means of a 
canonical transformation. In the relativistic case, this is possible due to a 2-dimensional 
remnant of the Pauli-Giirsey symmetry of massless fermions IM] and explains why 
two seemingly different large N field theory models give identical results [SSI- 
has to do is redefine particles into antiparticles for left-handed fermions only. If one 
works at non-zero chemical potential, a baryonic chemical potential /i in the GN model 
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corresponds to an "axial" chemical potential /is in the dual BCS-type model. The phase 
diagram which we have discussed in section 5 (for uiq = 0) is equivalent to the phase 
diagram of a theory with Lagrangian 

n 2 



N 2 



n=l 



N 



E 

Ln=l 



, (n)t , (n)t I (n.) I (n) 
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provided we reinterpret the chemical potential fi = fiR + fiL axial chemical potential 
= f'R — f-L- The kink-antikink phase of the GN model can then be identified with 
the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase [Sni EZ] of the dual model. Such 
inhomogeneous superconductors have recently attracted considerable attention in the 
context of QCD (for a review article, see reference [5Hj). 

Let us now turn to non-relativistic condensed matter physics and demonstrate 
that the GN model can also be related to quasi-one-dimensional superconductors in 
nature. In 1981, Mertsching and Fischbeck addressed the quasi-one-dimensional Peierls- 
Frohlich model with a nearly half- filled band, an electron-phonon system jHS]- This is 
the same basic model as the continuum model for degenerate polymers jH], but here 
the phase diagram at finite temperature was considered, notably the transition between 
commensurate-incommensurate charge density waves. Comparing with our results at 
mo = 0, we find a mathematical one-to-one correspondence between this system and the 
GN model now extended to finite temperature. The authors of reference jSH] have also 
found the analytic solution to the mean field equation, guided by the Landau expansion 
around the triple point (which is called Leung point [HDl in this context). 

In a subsequent paper, Machida and Nakanishi [HI] used the phase diagram 
of reference jHSj in a different physics context: They studied the interplay of 
superconductivity and ferromagnetism in ErRh4B4 (Erbium- Rhodium-Boride). They 
managed to reduce this problem mathematically to the Peierls-Frohlich model. For 
real order parameter, their results are again fully equivalent to ours for the GN model, 
except that now one has to use another dictionary: The Dirac equation corresponds 
to the Bogoliubov-deGennes (BdG) equation, flavor (A^ = 2) to spin (which exists in 
a quasi-one-dimensional world), chemical potential to magnetic field, baryon density to 
spin polarization. Our three phases (massive, crystal and massless) correspond to their 
BCS, "sn" and normal phases, respectively. Not only the phase boundaries, but all 
observables can be identified if one keeps in mind the above mentioned dictionary. Let 
us try to understand in more detail the reasons behind this remarkable correspondence. 
The BdG Hamiltonian involves four fermion fields, namely right (ips) and left (0^) 
moving electrons with spin up (s = -|-) and spin down (s = — ). This Hamiltonian can 
be mapped onto the Peierls-Frohlich Hamiltonian by the canonical transformation 



V 0- / 



V <Pl J 
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i.e., particle hole conjugation for spin-down fields followed by a spin-flip of left-moving 
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fields. Under this transformation, spin density goes over into ordinary fermion density, 
and one can understand all other relationships as well. 

Finally, we would like to mention the more recent work of Buzdin and Kachkachi 
[02] • They derive the Ginzburg-Landau theory for nonuniform LOFF superconductors 
near the tricritical point in the (T, if)-phase diagram in one, two and three dimensions. 
If we take their result for one dimension and specialize it to a real order parameter, we 
find perfect agreement between our equation and their equation (3) in appropriate 
units. Once again we have to identify their magnetic field Tio with our chemical potential 
/i for the reasons discussed above. 

12. Concluding remarks 

At the time of writing our previous review article on the thermodynamics of two- 
dimensional quantum field theories [HI, some progress had been made on fermionic 
models with continuous chiral symmetry like the 't Hooft or the NJL2 model. The fact 
that the existence of baryons implied a crystal structure of baryonic matter had been 
understood and the physical interpretation in terms of a gap at the Fermi surface — the 
"rediscovery of the Peierls effect in relativistic QFT" — had been clearly stated. Towards 
the end of this article we identified areas where future work was needed, mentioning in 
particular the GN model with discrete chiral symmetry. We hope to have shown in 
the present work that this problem has been solved in the meantime, including the 
generalization to the massive model with explicit symmetry breaking. As a result, the 
situation is now reversed: Today we know a lot more about the thermodynamics of the 
discrete chiral GN model than about the NJL2 model, particularly in the presence of a 
bare mass term. For a long time the phase diagram of the NJL2 model was thought to 
be identical to that of the GN model, but this can now be ruled out due to their different 
baryon structure. It will be interesting to see whether one can make further progress 
on the phase diagram of the massive NJL2 model where topology is expected to play a 
crucial role, analogous to the Skyrme model and Skyrme crystal in 3+1 dimensions. 

In reference we also pointed out that this type of soluble QFT models still had the 
potential to surprise us after so many years of studies. If anything, this impression has 
only been reinforced since then. The simple Lagrangian (P) has generated a much richer 
phase diagram than previously thought, together with some beautiful mathematics. 
Even more surprising for us was perhaps the discovery that a relativistic QFT "toy 
model" has such a close relationship to various quasi-one-dimensional condensed matter 
systems. We have tried to cover these less familiar aspects of the GN model in the present 
work as well. They represent a good example for a fruitful exchange between condensed 
matter and particle physics. Originally, condensed matter physics, in particular polymer 
physics, could profit of the particle physicists know-how on soliton-like baryons in the 
GN model to understand for instance electrical conductivity properties. Subsequently 
the polymer continuum models were vigorously developed, notably in the direction of 
crystal solutions where the expertise has traditionally been residing in condensed matter 



Updating the Gross-Neveu phase diagram 



32 



physics. During our work, we were able to take advantage of the progress in polymer 
physics to settle some unresolved issues in QFT, while at the same time generalizing 
the bipolaron lattice to finite temperature. 

This whole process was not as straightforward as it may sound. It was not that easy 
for us to find the relevant information in the unfamiliar literature where it was often 
hidden under a lot of material of less interest to us. Thus we could not avoid doing 
many things "the hard way" . For example the calculation of the phase diagram in the 
chiral limit, first numerically and then analytically, was performed by us independently 
of already existing pertinent results in the "parallel world" of condensed matter physics. 
In any case, our study of the massive model has profited immensely from this interplay, 
and it was better to realize and exploit the relationship between these two different 
branches of physics late than never. 

The author wishes to thank the guest editors of this special issue of J. Phys. A, in 
particular Gerald Dunne, for the invitation to write this contribution. A most pleasant 
and productive collaboration with Oliver Schnetz and Konrad Urlichs on the work 
reported here is gratefully acknowledged. I also thank both of my collaborators for 
a critical reading of the manuscript. 
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